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The finite sums of powers of cosecs occur in numerous situations, both 
physical and mathematical, examples being the Casimir effect, Renyi 
entropy, Verlinde’s formula and Dedekind sums. I here present some 
further discussion which consists mainly of a reprise of early work by 
H.M.Jeffery in 1862-64 which has fallen by the wayside and whose 
results are being reproduced up to the present day. The motivation 
is partly historical justice and partly that, because of the continuing 
appearance of the sums, his particular methods deserve re-exposure. 
For example, simple trigonometric generating functions are found and 
these have a field theoretic. Green function significance and I make a 
few comments in the topic of Renyi entropies. 
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1. Introduction 


The finite summations of powers of cosecs occur in a number of different ar¬ 
eas which, mathematically, are associated with image constructions, in one way or 
another, and involve the regular subdivisions or discretisations of the circle, a very 
old topic. 

Relatively recent papers that detail some of the history behind these, and like 
trigonometric summations, are Berndt and Yeap, [1], and Cvijovic and Srivastava, 
[2,3] and I refer to them for motivation, both physical and mathematical. I will 
not repeat the references in these useful works, unless they are directly relevant, 
but I might add a few more. These summations continue to appear in the physics 
literature. 

My intention here to give some historical and calculational details which I hope 
will be interesting and/or useful. I have elaborated the algebra because the main 
reference I employ (Jeffery, [4]) is, perhaps, a little obscure and certainly unknown. 
Also the often quite simple methods have application today. I have also included an 
example that is not in [4] as a slight novelty as it combines bosonic and fermionic 
elements. 

The next section introduces the sums in their general form and a (known) 
computable answer. Section 3 starts with Jeffery’s generating function for the 
simplest, untwisted sum for which an explicit expression is thence found. 


2. The summations 


The summations in question are the modulated cosec sums. 


t \ f 2 t:wI\ 2i/( 

C^,[n,w) = cos ^- j cosec ^ 


— ) , z/ e Z, 
n 


( 1 ) 


with twisting rc G Z, In fact here I will deal mostly with the untwisted case w = 0 
The sums were evaluated in terms of generalised Bernoulli polynomials by a 
contour method, [5] and [6], 


')2i' 


{n,w) = {w + iy\n,l), 


( 2 ) 


which can be expanded as a polynomial in n. 

This could be taken as the final answer, and it is, but it is a little awkward 
to expand hence I now pass to another expression, by a different approach, which 
yields a form amenable to hand calculation. 
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3. Jeffery’s generating fnnction 


In fairly recent times generating fnnctions for the (nntwisted) snmmations have 
been given by Fisher, [7], and Zagier, [8], (and, more generally, for the twisted ones 
in [6]) bnt the main point I wish to make in the present paper is that a mnch earlier 
version was derived by Jeffery, [4], in his discnssion of certain classes of integrals 
associated with the derivatives of the Gamma fnnction. 

Amongst other things, Jeffery derives the following resnlt in §36, ^ 

n —1 / \ -1 
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TT / TT , . TT , A f 

— > cot — [I — X) — cot —[l + X)]= / 

2n \ n n J Jo 

(y-^ - y‘) = f dy 
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>0 


=2 1 - 




JL 

1 - 
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1-r 

1-2/ 


y 


\y-^-y^) 


—xju _ ^.xln 


n(l - y) 




TT TTX 

= — cot- TT cot TTX , 

n n 


(3) 


where n can be even or odd. 

The third line resnlts from expanding y^ and then nsing, 

Jo _n-/0 




(-1)'^9!C(9 + 1) • 


I amplify Jeffery’s derivation of the hrst line of the identity, (3), which he does 
not give in detail. I will assnme that any formnla in (9) is known, althongh Jeffery 
derives many of them anew. Snbstitntion of the integral, 

into the standard resnlt. 


—TT cot TIZ — — '0(1 — z) 


gives 


—TT cot TIZ 



-t^-^ +t-^ 
1-t 


(4) 


^ I enlarge on his algebra. Be aware that there are a number of bad misprints in this tersely 
written, disjointed, but very interesting paper. 


3 



Therefore, z = (/ ± a:)/n, 


TT, ^ ^/7 \ 

TT cot —il — X) — TT cot —il + X) 
n n 

rl _j_ ^-{l+x)/rL _|_ ^(l-x)jn—\ _ ^-{l—x)lr 

= dt - 

Jo 1-t 

= n t dy 

Jo 


1-y^ 


where t = y 

Hence, writing ont the snm at length. 


n-T / ^ ^ \ pi 

TT } I cot — {I — x) — cot —{l + x)]=2n dy 

n J Jo 


l + y + ... + y 
1-y^ 


n—2 




which is the hrst line in (3). 

Eqnation (3) constitntes a generating fnnction and is nsed as snch by Jeffery 
to evalnate the snms 0 ) by hnding the (odd) derivatives with respect to x at 

0 . 

One might as well take the hrst derivative at once to give. 


1 /7r\ / 2'^ n \ 2'^ n \ 

“ ~ / cosec —(/ —xj + cosec —[l+x) 

2 \n J \ n 


n 


= 2(1-;^) C(2) + 2.3 (l - ^ j C(4) + 2.5 (l - ^ ) C(6) + ... 

2 2 2 '^^ 

= 1 TT cosec TTX - 5 - cosec — 

n 


(5) 


Setting X to zero yields the easiest snm. 


<^ 1 (^, 0 ) =- 5 ^ cosec f-j =—^C(2) = —^ 
1=1 ^ ' 


Jeffery’s paper is the earliest to which I can trace this resnlt. 

Eqnation (5) is nsed in the following way. The snms ( 72 ( 77 -, 0) and C^{n^Q) are 
given explicitly in [4] as farther examples and the constrnction of the nnmerical 
coefficients exhibited there indicates that [4] has employed, withont comment, the 
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recursion formula, 


cosec y = 


1 


(z/- l)(z/-2) 


[D‘1-\-{u — 2)“^^ cosec“^{y + z) 


, z/>2, (6) 


z=0 


for even z/, and = djdz, which must have been common knowledge, probably 
going back to Euler. Much later references are Ely, [13], and Saalschiitz, [14], 

The complete iteration of this yields. 


cosec = 


(2 p-l)! 


{Dl + (2p - 2)2) ... {Dl + 22) cosec 2 (y + z) 


z=0 


p-1 


(2P-1)1^0 


^Uf Df cosec2(y+ z) 


(7) 


z=0 


SO that equation (5) can now be employed to give the required even derivatives at 
zero. 


n—1 

DI^ ^ cosec 2 
Z=1 




2(2i + 1)! 



C(2i + 2) 

yj-2z-|-2 


Therefore, hnally, setting y = nl/n in (7) and summing. 



22p—2 

(2 p-l)! 


p-1 

E 


j=0 


{2i + 1) 

22i 


Vf 




C(2i + 2) 

7I-2i+2 ’ 


( 8 ) 


which I will refer to as Jeffery’s form. The constants Uf equal the sum of the 
products of the squared hrst p — 1 integers (excluding zero) taken p — 1 — i at a 
time, but, for later numerical convenience, I have extracted a factor of a power of 
two to give the more usual constants, W[. 

A quite recent paper, [15], has again considered these sums and reaches the 
same expression as (8) by a similar method. It appears also in a recent work, [16]. 


^ It is interesting to note that this recursion is an early manifestation of Hadamard’s technique, 
of increasing the dimension of the manifold by two by applying an intertwining operator to a 
propagation quantity. Hadamard worked in flat space. An application to spheres was made by 
ourselves, [10], and in a relevant conical context in [11]. See also [12]. 
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4. The calculation of the coefficients 


The combinatorial expression for the coefficients, Uf, directly applied, is satis¬ 
factory for a numerical hand calculation of the lower summations (Jeffery stops at 
p = 4) but is arduous for high p. For this reason consider the product occurring 
in (7). It is clear that the coefficient of equals times the coefficient of 

up to a sign, in the product 


— {p — 1)^) ... “ l) — ^ 

where is a even central factorial, see e.g. Stefiensen, [17]. I have called this 
coefficient, Wf. 

The expansion of the central factorial in terms of the central factorial coeffi¬ 
cients (essentially just a Taylor expansion), 


^[ 2 p ]-2 


P jj2u Q[2p] 


is here expressed as central derivatives of nothing. These satisfy a recursion which 
allows machine calculation. This is, [17], 


]j2u q[2p+2] jj2u-2 q[2p] ^ jj2iy q[2p] 

“ {2u-2)l {2u)l 


with the starting values 

0[2p] 


( 2 )! 


(-ir-‘((fc-i)y; 


]j2p q[2p] 

(2p)! 


Hence the coefficients in (7) are, ^ 


W[ 




]j2i' q[2p] 


(9) 


Early tabulations can be found in Stefiensen, [17], and Thiele, [19] p.35 suf¬ 
ficient to reach p = 8. Later tabulations exist (e.g. in [20]). Table 1 in the more 
recent discussion of these cosec sums by Grabner and Prodinger, [16], is the same 
as the one in Thiele and equivalent to that in [15]. 

^ Jeffery, [18], gives a symbolic form for the combinatorial description but it is expressed in terms 
of forward differences and is more complicated. 
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I table some values taken from Thiele, 


p 1 2 3 

4 

5 

i 

1 1 4 

36 

576 

0 

1 5 

49 

820 

1 

1 

14 

273 

2 


1 

30 

3 



1 

4. 

Substitution of these numbers into 

the 

answer. 

(8), produces agreement with 


existing results and, of course, with (2). There is no need to write out any specihc 
cases. 


5. Jeffery’s twisted generating function 


Jeffery also gives expressions for the (fermionic) summations,^ 

" ’(-l)‘coT-Tosec“(- 

\ n 

/^i\ 

. 2m-\-l 


- = ^(-l)^cos ( ■ 

1=1 ^ 

= ^(—l)^cot ^ cosec(^— ) , 
1=1 ^ ^ 


V n 


( 10 ) 


which, again, he derives from a generating function, 

n—1 

—1)^"*“^ ( cosec 

n ' ' n 


n-i , X 

TT V—v . / I \ n \\ 

— > (~l) cosec — [I — X) — cosec —[l + x) = n cosec ttx -cosec — , (ll) 

2n \ n n J n n 


if n is odd. 

The proof proceeds as for the untwisted case and begins with the representation, 

pi ^ 


TT cosec TTo: = dt 

Jo 


I _ ^-x 

... =“+/ , 1 ^ 1^1 < 1 , ( 12 ) 

t + 1 X Jq t + 1 


which follows by differentiating Rummer’s integral, (re-proved by Jeffery), 


, TTx 

logtan— = dy - 

2 Jo (1+ 2 /) logy 


^ I give with elaborations the, sometimes elementary, algebra because the reference may not be 
available. The journal can actually be found in the Gottinger Digitalisierungscentrum. 
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Hence, 


TT cosec — (/=F 3:) = / dt 
n In 


if I lies between 0 and n 

Set t = y~^ then the integral leads to, 


1 f(lT^)/n-l _ ^-(lTx)/r 


, |a:| < 1, 


— ( cosec — {I — x) — cosec — (/ + x) ) = 


- / dy^ -^^^- (13) 

Jo 1 +r 

_ . 1 ^ (^y-l-l+n _ yl-l'j (^y-X _ ^ ^ 

~ Jo ^ i + r 

This allows the snmmation to be done. 

Now pnt in the valnes of / = 1,2, ...,n — 1 in (11) with the correct parity. The 
nnmerator of (13) contains (n is odd). 


n—2 ^n—3 


y^-^ + ...-y + l + l-y + y^ - ... + y^-^ 


so that (13) becomes. 


^ ( cosec ^(/ -x)- cosec ^(/ + x)] = 


= / dy 


= / dy 


= / dy 


{1-y + y^ - ... + y^ ^){y ^ - y^) 


l + y'^ 

_1 - 

l + y l+y^ ^ ’ 

y-X _ yX \ y-xju _ yXjj 

1 + y n 1 + y 


TT 'KX 

= TT cosec Tvx -cosec — 

n n 


after nsing the second representation in (12) for cosec. We have proved (11). 

The procednre is exactly as previonsly. The power series of the right-hand side 
of eqnation (11) is obtained by expanding y^ and nsing. 


i' - Ei‘ (-1)^ ‘ H-DN! 


oo „1 


{-ly q\y{q + l) , 



where rj is the Dirichlet ry-function, related to the Riemann ^-function by 77(5 +1) = 
(1 — 2“'3') C,{q + 1). (Note that the trigonometric closed form of the generating 
fnnction has not actnally been nsed.) 

This resnlts in, 


TT 

2 n 


71 1 / \ 00 / 1 \ 

/ cosec —{I — x) — cosec —{l + x) \ =2 ^ ( 1- — j , 

which I differentiate with respect to x to get, 

- f^ (—( cot —{I — x) cosec —{I — x) + cot —{I + x) cosec —{I + 

2 \n J \ n n n n J 


J = 1 


(15) 


For the left-hand side, note that 


cosec {y -\- z) = —{2m + 1 ) cosec {y + x) cot{y + z) (16) 
and the iteration for cosec cot is determined by that for cosec which is. 


cosec 


2m+l,, _ 


y = ~ • • • (^^ + 1 ) cosec {y ± z) 


(2m)! 


(17) 


z=0 


and can again be expanded in terms of central factorial nnmbers. I write ont the 
expansion in the form. 


cosec ^^+^y = 


1 


(2m)! 


Vl^Dl^cosec {y ± z) 


i=0 


(18) 


z=0 


Differentiating (18) with respect to y, nsing (16), setting y = 7r//n, z = ±7ra:/n 
and taking the snm. 


n — 1 


ttI 




cosec 


2m+l 


tt/ 


1=1 


E'".” 


(2m+ 1)! Y 


n 


TT 


n 

2i n — 1 


1=1 


(“1)^"*"^ f cot —{I ± x) cosec —{l±x) 

f ^ \ fl Ti 


fr=0 
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Substituting in the derivatives at a: = 0, which can be read off from (11), 
delivers the final result, for n odd. 


5^(-l)'+icot 


^ ’ i=0 

which should be compared with the bosonic sum, (8). One notices now the appear¬ 
ance of the Dirichlet function, typical for fermionic quantities. 

From the form of the full iteration, (17), the coefficients, are given by 
the sums of products of the squares of the odd natural numbers, and this is how 
Jeffery, [4] §41, computes them. Another way is to employ a recursion relation but 
the easiest option is to use Thiele, [19], p.36 who conveniently tabulates them as 
the positive integers they are. These are sufficient to reach m = 8. 

For the convenience of the reader I lift a few lower values from [19], 

m 0 1 2 3 4 i 

1 1 9 225 11025 0 

1 10 259 12916 1 

1 35 1974 2 

1 84 3 


The tables in Steffensen, [17], and in [20] can also be consulted. 

As an example, Jeffery writes out the specific case m = 2 and I here repeat it, 
putting in the numerical values of the coefficients to give an explicit polynomial. 


5^(-l)'+icot 


— cosec°— = — (n^ ~ 1) H-(^^ ~ 1) J-(^° ~ 1) ? 

n n 80^ ^720^ ^ 15120^ 


if n is odd. 


6. A mixed summation of alternating cosecs 


To ring the changes I compute, using Jeffery’s method, the alternating sum 


-C'^(n,n/2) = 


cosec 


, e z. 
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for even n, which is neither bosonic nor fermionic. The algebra is a mixtnre of that 
in sections 3 and 5 and starts with the representation (4) for the cotan, so that 


— ( cot —{I — x) — cot —{l + x) 

2,n \ r) n 


n—1 


1=1 


= / dy 


= / dy 




TT 

n 

n—2 


1 y 


l-yn 

n — 1 


(y-^-y^) 


+ 


1 + 2/ 1 -y 


= 2(^r7(2) + ^Ja:+ 2(77(4) + 

TT Tra: 

= TT cosec TTX -cot - . 

n n 


(y"' - y') 


C(4)\ 

77-"^ / 


dy 


y — ^ _ y—xjn _ yXjr 

l + y ^ 


77(1 - y) 


X + . . . 


( 21 ) 


Differentiating with respect to x, 

2 n—1 

V 77 


- f ^ (—( cosec^ —(/ — x) + cosec^ —(/ + x) 

2 \Tl J \ n n 


1=1 


77 


1 


= 2 77(2) + ^C(2) +2.3 


77^ 


‘(>7(4) + +2.5(^ 17 ( 6 ) + ^C( 6 )^ +“ + ... 

( 22 ) 


Again, the recnrsion for cosec ‘^^x is bronght in and gives, as before, for 77 even. 


- y(-l)^+^ cosec 2P — 

2 ^ 77 

= t ^ + 2) + C(27 + 2)) 


1 

^2i+2 ’ 


(23) 


which has both bosonic and fermionic aspects. The pnrely fermionic snmmation, 
(10), contains an extra ‘phase factor’ which is jnst the trace of the spin-1/2 rotation 
matrix throngh angle 27r/77 which is necessary becanse of the rotation of the zwei- 
beine, [21]. 

Formnla (23) is easily programmed. Nnmerical examples agree with the general 
expression, (2). Chn and Marini have compnted some examples, [22] p.l49, nsing 
a method somewhat less convenient than the one employed here, and I have fonnd 
agreement. 
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7. Comments, connections and conclusion 


As mentioned in the Introduction, the polynomials crop up in several seem¬ 
ingly disparate topics, examples being statistical mechanics, vector bundle theory, 
[6], Dedekind sums, interpolation, quantum held theory on multiply connected man¬ 
ifolds, such as lens spaces, or on orbifolds with conical structures, real or artihcial. 

In the latter category is the construction of the entanglement and Renyi en¬ 
tropies using the replica technique which introduces a conical singularity into a 
codimension 2 submanifold. This process is usually accomplished by using integer 
coverings but can equally well be done with images as is well known. In this connec¬ 
tion, Cardy, [23], has also recently recalculated the classic summations, Ci{n^ 0) and 
( 72 ( 77 -, 0). Furthermore in a related subject, Herzog and Nian, [24], have employed 
the recursion (6). 

In this regard it is interesting to note that the simple rewriting of the generating 
function (5),® 

n—1 

TT 

cosec ^ — ± x) 

(24) 

= 77^ cosec , 

gives the Green function in four dimensions as in [24] equn.(33), to a constant factor. 

More general summations (in higher dimensions) of the same form appear in 
the appendix to [6]. 

It should be remarked that, for the requirements of the present paper, the 
construction of generating functions such as (3) is somewhat of an afterthought. 
Only the, prior derived, power series are used. 

The relation between coverings and images can be illustrated by Lubbock’s 
summation formulae, [17], in interpolation theory as described in [26]. All I point 
out here is that the Lubbock approximation involves a series of polynomials which, 
in one variant, are shown, [26], by contours, to be given by, 

P 2 M = +(^ -1)/2 I 1) , 

in terms of generalised Bernoulli polynomials, or, 

P2A'n-) = ^^C^{n,{n-l)/2) , (25) 

® This is a by now standard identity, [25] p.211. 
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in terms of the fermionic snmmations/ Here, n wonld have to be an integer for the 
snm to make sense. After the snmmation has been effected, n conld be anything and 
in the context of Lnbbock’s formnla n is set eqnal to l/h with h integral providing 
the alternative form (see Steffensen, [17]) in terms of central factorials, 

{h-l)/2 . ..[2u] 

E ' 


^’2.(1/A) = 


1 


{2u)\ 


-(h-l)/2 


where the snm is over j in steps of one. 

Steffensen lists a few examples of the P. These check with (25) which is no 
snrprise in view of the image relation. 


(a + i I 1) = fe (a I 1.1). (26) 

s=0 

It might be helpfnl to spell things ont. The left-hand side is the pre-image snm 
giving the qnantity (here B) on the (mnltiply connected) circle of circu mf erence 
l/h in terms of that on a (bigger) covering ® circle of nnit circu mf erence consisting 
of h copies of the smaller circle in the form of circu mf erence intervals of length 1/h 
with interval ends. This can be referred to as a subdivision of the bigger circle. 

By contrast, an integer n-cover of the nnit circle means n copies of the nnrolled 
nnit circle joined end to end and the two bonndary points identihed giving a wrapped 
np circle of circu mf erence h. The left-hand side of the snm, (26), then does not 
make sense bnt the right-hand side does. 

When embedded in the plane, all circles have the same radins and give conical 
singnlarities with positive and negative dehcit angles for snbdivisions and integer 
coverings respectively. 

The specihc snmmations considered here all fall into the general twisted class, 

(1) , and if one is looking for jnst a nnmerical polynomial then the general expression, 

(2) , in terms of Bernonlli polynomials provides it. Jeffery’s forms, (8) and (19) 
(and (23) ) arrange these polynomials in a different, more explicit manner. This 
arrangement has advantages when considering the heat-kernel coefficients in the 
presence of conical singnlarities, [27]. 

Everything that has been given can be transcribed by replacing cot by tan and 
cosec by sec becanse sec and cosec obey the same recnrsion. (Some adjnstment of 
the snmmation is needed.) 


^ Another Lubbock variant gives polynomials, Q, which are the untwisted sums. 
® The term ‘covering’ is meant in the projection sense. 
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